Dipolar ordering in crystals of Mni 2 Ac 
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Ordering in realistic elongated box-shape crystals of the molecular magnet Mni2 Ac is investigated 
with the site-resolved mean-field approximation that does not assume a uniform ordering. It is shown 
that ferromagnetic ordering should not occur in crystals with the aspect ratio up to 12. Instead, 
for the aspect ratio about 6, the inner and outer regions of the crystal order in different directions, 
uniformly along its length. Finding ordering temperature by extrapolating the inverse susceptibility 
curve does not provide a correct To- 
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PACS numbers: 75.50.Xx, 75.30.Kz, 75.50.Lk, 64.70.Tg 



Magnetic ordering in molecular magnets such as Mni2 
Acji attracts attention of researchers because the order- 
ing dynamics, as well as domain-wall motion^ might 
be intimately related to the spectacular phenomenon of 
resonance spin tunneling 3 under the energy barrier cre- 
ated by the strong uniaxial anisotropy 4 - acting on the 
effective spin S = 10 of the magnetic molecule. The 
anisotropy barrier is responsible for the bistability of 
magnetic molecules, so that at temperatures below 10 K 
the latter can be considered as Ising spins 1/2 pointing 
up or down. Mni2 Ac crystallizes into a body-centered 
tetragonal lattice with the periods a — b — 17.319 A and 
c = 12.388 A, c being the direction of the magnetic easy 
axis z. The unit-cell volume is vq = abc = 3716 A 3 and 
there are two molecules per unit cell, v = 2. 

Since the magnetic core of magnetic molecules is sur- 
rounded by organic lygands, there is practically no ex- 
change interaction between the molecules. Thus the sole 
source of magnetic ordering is the dipole-dipole interac- 
tion (DDI) that leads to ordering below 1 K in Fes (Refs. 
HI!), Mn i2 Ac (Ref. @), and other compounds.^— The 
type of dipolar magnetic ordering in molecular magnets 
is a subtle question. Since a column of spins up directed 
along the c axis produces a magnetic field up on its own 
spins that largely exceeds the field it produces on other 
spins^ one concludes that the ground state includes fer- 
romagnetically ordered columns along the c axis. The 
ferromagnetic columns can order ferro- or antiferromag- 
netically with respect to each other. 

It was shown that the ground state of ellipsoids of rev- 
olution depends on their aspect ratio^ For ellipsoids very 
prolate along the c axis, the ground state was found to 
be ferromagnetic, whereas in other cases there are fer- 
romagnetic columns or planes ordered antiferromagneti- 
cally with respect to each other. The energies of differ- 
ently ordered states are rather close, so that practically 
the system will order in a spin-glass state upon fast cool- 
ing. 

Evidence for a ferromagnetic ordering in Mni2 Ac was 
obtained in Ref. by neutron scattering experiments. 
However, the shape of the crystal was not specified. Usu- 
ally crystals of molecular magnets are box-shaped and 
elongated along the c axis but it is a question if they are 
elongated enough to make ferromagnetic ordering prevail. 
Recent measurements of the ferromagnetic susceptibility 



X above T c in Mni2 Ac and its modifications 1 ^ were done 
on box-shape single crystals with the aspect ratio L z / L x 
of about 5-6. In the absense of transverse field, one can- 
not approach T c because the relaxation time becomes 
too long below the blocking temperature 3 K due to the 
anisotropy barrier. In this case T c ~ 0.9 K was obtained 
by the linear extrapolation of the x _1 (T) curve. 

Ap plying a strong transverse field B± up to 6 T (Refs. 
0G3) increases the relaxation rate via barrier lowering 
and spin tunneling, so that ordering can be achieved dur- 
ing a realistic timeZ and the susceptibility measurements 
can be extended to lower temperatures- 12 However, the 
transverse field tends to suppress magnetic ordering and 
lower T c via the three main effects: (i) spin canting that 
reduces the effective magnetic moment along the z axis, 
(ii) tunneling hybridization of the |±5) states that acts 
as the transverse field in the Ising model that lowers 
T c and leads to quantum criticality, and (iii) random 
longitudinal fields resulting from the strong transverse 
field and random tilts of the molecular easy axes, further 
suppressing ordering. These effects have been recently 
discusse d 13 ' 14 in connection with experiments 1 ^ on the 
basis of the mean-field approximation (MFA) that should 
work well for long-range interactions such as DDI. 

The aim of this work is to investigate the possibility 
of ferromagnetic ordering in its competition with other 
ordering types in crystals of Mni2 Ac of the realistic 
box shape. While calculations of dipolar fields in crys- 
tals of molecular magnets were done for ellipsoids of 
revolutio n 2 i 6 i 13 that are impossible to grow, no theoreti- 
cal work has been done yet on box-shape crystals. It is 
not obvious that long rods behave similarly to long ellip- 
soids. It was shown that the dipolar field in the middle of 
the end faces of a uniformly magnetized cylinder of Mni2 
Ac is opposite to the magnetization (the top of left col- 
umn of page 4 of Ref. [2) . This should lead to spin flips at 
the ends of the cylinder with a subsequent proliferation 
into its body. The same can be expected for elongated 
boxes. Since the dipolar field in crystals of other than el- 
lipsoidal shape is non-uniform, the MFA equations take 
the form of a large system of equations for all magnetic 
molecules considered separately. 

The model includes pseudospin variables cr, = ±1 for 
molecules at each lattice site i of a boby-centered tetrago- 
nal lattice. The magnetic moment of a molecule is Sgfi B 
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with g = 2. The dipolar field on magnetic molecule i is 
the sum over positions of all other molecules j 
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Here D zz is the reduced dipolar field and 
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Inside a uniformly magnetized ellipsoid, a z = const, the 
dipolar field is uniform and one has D zz = D zz a z , where 
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is the number of magnetic molecules per unit cell (y = 2 
for Mni2 Ac) and — 0, 1/3, and 1 for a cylinder, 
sphere, and disc, respectively. The reduced dipolar field 

in a sphere Z?iz ph ' ) depends on the lattice structure. For 
Mni2 Ac direct lattice summation yields D z s z ph ^ = 2.155 



that results in D 



(cyl) 



10.53 for a cylinder^ Then Eq. 



52.6 mT in an elon- 



(JTJ yields the dipolar field B 
gated sample that was also obtained experimentally. 15 
The ground-state energy in the above uniform states is 
given by 



E a = -(l/2)D zz E L 



(Sg» B f/v , (4) 



where Ed is the dipolar energy, EdI^b — 0.0671 K for 
Mni2 Ac. The mean-field Curie temperature is given by 2 - 



T c = E D D„/ki 
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that results in T c — 0.707 K. 

States with ferromagnetically ordered planes alternat- 
ing in the a or b directions in each sublattice of Mni2 Ac 
have D zz = 9.480, independently of the sample shape. 2 
The state with alternating chains in each sublattice, 
directed along the c direction has a very close value 
D zz — 9.475. For the two-sublattice antiferromagnetic or- 
dering one obtains D zz — 8.102. Thus, in a strongly pro- 
late ellipsoid of Mni2 Ac ferromagnetic ordering is pre- 
ferred. It is interesting to estimate how strongly prolate 
the ellipsoid has to be for this to be the case. Equating 
D zz of Eq. ([3]) to D zz = 9.480 for the alternating-planes 
structure, one obtains the maximal demagnetizing factor 
rS^ — 0.0419. Using the formula for prolate ellipsoids 
of revolution, one obtains that the minimal shape aspect 
ratio R z /R x = 6.13 is required for ferromagnetic order- 
ing. 

The longest crystal used in the experiments of Ref. [l2| 
has dimensions 0.4 x 0.4 x 2.4 mm 3 and thus the aspect 
ratio 6 that would be still slightly insufficient for a crys- 
tal of ellipsoidal shape to order ferromagnetically. We 
will see below that even much longer box-shape Mni2 Ac 
crystals do not order ferromagnetically. The reason for 
this is the above mentioned instability of the ferromag- 
netic ordering at the ends driven by the negative value 
of the dipolar field, D zz = -2.03,2 



The system of Curie- Weiss equations for a crystal of 
molecular magnet in a uniform external field B z has the 
form 
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where h z = Sgfi B B z /k B - The linearized Curie- Weiss 
equations above T c can be cast into the matrix form 



(TI - V) • (tr,) = h z l. 



(7) 



Here V is the DDI matrix, {V} 4 . = {E D /k B )4>ij, I is a 
unit matrix, {1}^ = 6%j, and I is a unit vector, = 1. 
In fact, position of a molecule in the lattice is defined by 
three numbers i a , %, and i c corresponding to the three di- 
rections a, b, and c, plus the sublattice index. The indices 
i,j are compound indices running from 1 to the number 
of sites in the lattice N. In terms of i,j the matrix V is 
non-symmetric, so that care should be taken by distin- 
guishing between its right and left eigenvectors. One can 
seek the solution in the form (cr z ) = ^Z^C^A.^, where 

are right eigenvectors, V • = T M A^, and T M are 
eigenvalues, fi = 1, . . . , N. The right and left eigenvectors 
satisfy the orthonormality condition A^, • A^ = 5^^, 
i.e., the matrix of left eigenvectors is the inverse of that 
of the right eigenvectors. Using the orthonormality, from 
Eq. one obtains C M = (A^ • I) hj (T — T^) . For the 
susceptibility of the whole crystal x = (1/^0 ((°"z) • I) /h 
one finally obtains 
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The ordering temperature Tc can be identified with the 
maximal eigenvalue T M , this is the first time as Eq. ((5J) 
produces infinity as temperature is lowered. 

Numerical calculations have been done using Wol- 
fram Mathematica for box-shape crystals of dimensions 
L a = Lb <C L c in units of the lattice spacing a. The 
numbers L a = and L c have been taken even. The 
indices i a and i& specifying positions of molecules within 
the a, b plane run in the range -i a ,6, m ax < ia.b < V6,max 
for sublattice 1 and in the range — « a .6,max + 1/2 < 
i a .b < Vfc,max — 1/2 for the body-centered sublattice 2, 
where i a ,b,max — L a /2. For the crystallographic direc- 
tion c, the index ranges are — i c ,max < «c < Vmax and 
-i c ,max + 1/2 < i c < «c,max - 1/2, respectively, where 
*c,max = Round \L a j (2^)] and n — c/a — 0.7153 for 
Mn 12 Ac. The total number of molecules in the crystals 
studied thus is given by = (2i Qiinax + 1) (2i b , max + 1) 
(2i c ,max + 1) + (2vmax) (2ib,max) (2i c , ma x) that reached 
N = 31061 for the biggest crystal with L a = Lb = 10 and 
L c = 100. However, the actial number of lattice sites in 
the calculations was reduces by a factor of approximately 
8 by using symmetry that resulted in the gain of approxi- 
mately 8 2 = 64 in the computer time and memory usage. 
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FIG. 1: Inverse susceptibility of elongated box-shape crys- 
tals of Mni2 Ac vs temperature. While the x _1 (T) curves 
extrapolate to the would-be ferromagnetic ordering tempera- 
tures, the real ordering occurs at a higher temperature and is 
non- ferromagnetic . 



FIG. 2: Magnetization of elongated box-shape crystals of 
Mni2 Ac of different dimensions, developing upon lowering 
temperature at a slow constant rate. Static inverse suscepti- 
bility for L a — Lb — 8, L c — 50 from the preceding figure is 
also shown. 



Numerical calculations show that the uniform-field 
coupling coefficient (A^ • I) (A^ • I) reaches a large 
value for some /i = fJ-f erro , while for other /i values 
(A^ • I) (A^ • I) are much smaller. The ferromagnetic 

eigenvectors A^ R = A.f^ a have all their elements of 
the same sign, while other eigenvectors have elements of 
different signs and do not project well on the uniform 
field. Thus not too close to Tc Eq. © is dominated by 
the single ferromagnetic term, 

g 1 ( A fcrro • jj ( A fcrro [ jj /g-, 

N T-T ferro 1 ' 

and % _1 (T) cc T — Tf crro is a straight line that extrapo- 
lates to Tf erro as an apparent transition temperature. In 
all cases for which numerical calculations have been per- 
formed, it turns out that Tf crro < Tc and thus Tc cor- 
responds to other types of ordering than ferromagnetic. 
At T approaches T c , the term with T M = T c in Eq. © 
becomes dominant and the curve x _1 (T) drops suddenly 
to zero. The ordering eigenvector corresponding to the 
ordering at Tq was shown to be non-ferromagnetic. Of 
course, the behavior shown in Fig. [1] for several different 
crystals cannot be seen and actual Tc cannot be found, 
if only the high-temperature susceptibility data is avail- 
able, as in Ref. HH In fact, there are many different 
ordering eigenvalues in the region between Tf crro and 
T c 

It appears that with increasing the crystal size cou- 
pling of the uniform field to non-ferromagnetic eigenvec- 
tors decreases, so that the drop of at Tc be- 
comes sharper. This makes using uniform susceptibility 
to detect Tc in the case of a non-ferromagnetic ordering 
questionable. Of course, theoretically one can couple to 
the non-ferromagmetic ordering modes by a non-uniform 
magnetic field (such as the staggered field in the case of 
antiferromagnetism) but practically it is difficult to real- 
ize. 

One can see from Fig. Q] that Tf orro essentially depends 
on the aspect ratio, as it should be, and is in accord with 



Eq. j3|) . To the contrary, Tc does not strongly depend on 
the shape since it corresponds to a non-ferromagnetic or- 
dering. Its moderate increase with the crystal size shows 
that the crystal sizes in the numerical calculations are 
still somewhat small to perfectly reproduce the behavior 
of macroscopic crystals. With increasing elongation, fer- 
romagnetic ordering becomes more competitive but still 
falls behind other orderings up to the aspect ratios of 
about 12 where ellipsoids would already order ferromag- 
netically. It is difficult to increase the elongation while 
keeping the transverse size large enough in the calcula- 
tions because the number of molecules becomes too large. 

Note that in the case of the standard antiferromag- 
netism the ferromagnetic eigenvalue Tforro is negative. 
The ferromagnetic state in the antiferromagnet is abso- 
lutely unstable, as the molecular field is opposite to the 
spins. In our problem of dipolar ordering, there are many 
positive and many negative eigenvalues T M . Negative 
correspond to absolutely unstable states, whereas Ij, > 
correspond to local minima of energy with dipolar fields 
parallel to spins. 

Competition of many local energy minima in our model 
of ordering in Mni2 Ac makes it impossible to describe 
the ordered state below Tc by solving Eq. ^ directly. 
Instead of the global minimum of the free energy, the 
solver finds local minima, local maxima, or saddle points. 
A more reliable method is to solve the system of relax- 
ational equations 

d (a lz ) /dt = -T [(a tz ) - tanh (. . .)] , (10) 

where T is the relaxation rate and tanh (...) is the same 
as in Eq. ©. As the purpose of this work was to study 
ordering rather than the exact dynamics, T was set to 
an arbitrary constant in numerical calculations. It was 
found that Eq. (fTOf leads to freezing into a spin-glass 
state as the result of relaxation out of a random initial 
state at low temperatures. To obtain the ordering type 
that is mostly close to the lowest free energy state at any 
temperature, one can solve Eq. (|10p with temperature T 
slowly changing in time from some T max > Tc to nearly 
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FIG. 3: Ordering in elongated box shape Mni2 Ac crystals 
(L a = Lb = 8, La = 50) at T = 0, obtained by slow lowering 
temperature from T max > Tc and shown in the cross-section 
through the middle of the crystal, a) Sublattice 1; b) Sub- 
lattice 2. Magnetization in both sublattices is qualitatively 
similar. Central regions are magnetized opposite to the small 
external field while periferal regions are magnetized along it. 
Magnetization at T — is uniform along the c direction. 



zero. Just to initiate ordering, one can set h z to a very 
small value. 

Fig. [2] shows the results of these calculations for two 
different crystal sizes. The (a z ) curves accurately repro- 
duce the equilibrium magnetization everywhere except 
for the critical region (0.6 K < T < 0.7 K) where critical 
slowing down requires a slower temperature sweep that is 
difficult to implement in the numeric routine. The aver- 
age magnetization at T — > in both cases is significantly 
smaller than 1 because of a non-ferromagnetic type of or- 
dering. A nonzero average magnetization in the ordered 
state seems to be a finite-size effect related to the incom- 
plete compensation of molecules with spins up and down 
in small crystals. Ordering attained at T — > is shown in 
Fig. [3] for cross-sections through the middle of the crystal 
length. In fact, the magnetization at T — > is uniform 
along the c axis. The inner and outer regions of the crys- 
tal order in different directions, in a similar way in both 
sublattices. 

Concluding, mean-field calculations do not support fer- 
romagnetic ordering in elongated box-shape crystals of 
Mni2 Ac that are currently under investigation. 12 How 
long must be box-shape crystals to order ferromagneti- 
cally remains an open question. In any case, ferromag- 
netic ordering should be extremely difficult to observe 
in experiments because of a strong competition from nu- 
merous other types of ordering. 
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